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ABSTRACT 

The detailed distribution of mass within clusters of galaxies can be used to probe the nature of dark 
matter. We show that constraints on the extent of the mass distribution around galaxies in the rich 
cluster A 2218 obtained from combining strong and weak lensing observations are consistent with the 
predictions which assume that the dominant mass component (dark matter) in these halos is coUisionless. 
A strongly interacting (fluid-like) dark matter is ruled-out at a confidence level of more than ha. 

Subject headings: gravitational lensing, galaxies: fundamental parameters, halos, methods: numerical 



1. INTRODUCTION 

The detection of gravitational lensing of distant galax- 
ies by foreground galaxies has been used in recent years 
to set constraints on the masses and sizes of the lensing 
galaxy halos. Current studies suggest the presence of halos 
extending beyond 100 kpc (Brainerd, Blandford & Small 
1996; Griffiths et al. 1996; DeU'Antonio & Tyson 1996; Fis- 
cher et al. 1999; McKay et al. 2001; Wilson et al. 2001). 
The same technique can also be applied to galaxies in the 
dense environments of massive X-ray clusters. Analysis 
of galaxy-galaxy lensing in the cores of rich clusters im- 
plies a smaller spatial extent for the dark matter halos 
associated with morphologically-classified, early-type clus- 
ter members than that seen in the field (Natarajan et al. 
1998). Most recently, Natarajan, Kneib, & Small (2002) 
[NKS02 hereafter] have found compelling evidence for the 
tidal truncation of dark matter sub-halos in clusters based 
on Hubble Space Telescope imaging of lensing clusters at 
intermediate redshifts. 

The detailed mass distribution within lensing X-ray 
clusters, and in particular the fraction of the total clus- 
ter mass that is associated with individual galaxies - has 
important implications for the frequency and nature of 
galaxy interactions. As we show in this Letter, the char- 
acteristic halo sizes of cluster galaxies that survive tidal 
deformation and stripping in their dense environments of- 
fer tantalizing clues as to the nature of dark matter itself. 
While there is incontrovertible evidence for the existence 
of significant amounts of dark matter in the Universe, its 
nature remains one of the challenging and unsolved prob- 
lems in cosmology. Following the original suggestion by 
Furlanetto & Loeb (2002), we demonstrate here that con- 
straints on the properties of the dark matter, whether it is 
coUisionless or collisional (fluid-like), can be derived based 
on the truncation radii of galaxy halos in clusters. 

If the dark matter is coUisionless, then the sizes of galaxy 
halos in clusters are primarily shaped by the dynamical ef- 
fects of tidal truncation and collisional stripping (Ghigna 
et al. 1998). The dominant process of tidal truncation acts 



on the short orbital time-scale due to the tidal field of the 
cluster as a whole (Taylor & Babul 2001), while collisional 
stripping results from binary interactions among individ- 
ual galaxies (Binney & Tremaine 1987). The global tidal 
field of the cluster truncates the dark matter halo of each 
galaxy at a radius inside of which the mean mass density 
of the galaxy is roughly equal to the mean interior density 
of the cluster. However, if the dark matter is fluid-like 
then galaxy halos would be further stripped by ram pres- 
sure down to a radius which is significantly smaller than 
the tidal radius (roughly by the ratio between the velocity 
dispersion inside the galaxy halo and the velocity of the 
galaxy through the cluster; see Furlanetto & Loeb 2002 
for further details). In this work, we use the truncation 
radii of a sample of cluster galaxies, inferred observation- 
ally from the galaxy-galaxy lensing in the cluster A 2218 
(NKS02), to decide whether the dark matter is coUisionless 
or fluid-like. 

2. USING GALAXY-GALAXY LENSING TO MEASURE 
TRUNCATION RADII OF CLUSTER GALAXIES 

We model the background cluster as a superposition 
of a smooth large-scale potential 250" in extent) and 
small-scale clumps that are associated with bright early- 
type cluster galaxies (see Natarajan & Kneib 1997 and 
Natarajan et al. 1998, for more details). To quantify the 
lensing distortion induced by the global potential, both 
the smooth component and individual galaxy-scale ha- 
los are modeled self-similarly using a surface density pro- 
file, S(i?), which is a linear superposition of two pseudo- 
isothermal elliptical components (see the PIEMD models 
derived by Kassiola & Kovner 1993), 



S(i?) 



1 



1 



1 - roln \ ^rl + i?2 ^^fTW 



(1) 



with a core radius tq and a truncation radius r* 3> tq. 
The free parameters of this profile are chosen for both the 
smooth component and the clumps so as to obtain the ap- 
propriate mass distributions on the relevant scales. The 
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projected radius R is a function of the sky coordinates x 
and y and the elhpticity e (see §2.2 of Natarajan & Kncib 
1997). One of the attractive features of this model is that 
the total mass is finite (oc Eorort). With the additional 
assumption that Ught traces mass, galaxy halos in clusters 
are characterized by the following scaling laws: 

'^ = "*(£)'''"° = ^°*(£)^ = (I:) -^'^ 

These imply the following scaling for the ratio 

a-l/2 
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The total mass Mtot(oo) then scales with the luminosity 
as 



-Mtot(oo) = 27rI]orort = -^n = ^ct* n* I — 
and the mass-to-light ratio T is given by 
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For a = 0.5 the model has constant T for all galaxies (al- 
though T may be a fimction of radius inside each galaxy) . 
If a > 0.5 {a < 0.5) then brighter galaxies have a larger 
(smaller) T than fainter ones. The parameters that char- 
acterize both the global smooth component and the clumps 
are optimized using the observational data (i.e. the po- 
sitions, magnitudes, geometry of strong lensing features, 
and the smoothed shear field) as constraints. We find that 
a = 0.5, is the favored value for the best-fit mass models 
and in fact, a < 0.5 is ruled out at a very high confidence. 

Our goal is to optimally partition the total cluster mass 
between a smooth component and the clumps based on 
the observational data. A maximum-likelihood method is 
used to obtain significance bounds on fiducial parameters 
that characterize a typical halo in the cluster. The 
likelihood function of the estimated probability distribu- 
tion of the source ellipticities is maximized for a set of 
model parameters, given a functional form of the intrinsic 
ellipticity distribution measured for the faint background 
galaxies. For each faint galaxy j with measured shape 
Tobs) the intrinsic shape rg^ can be estimated in the weak 
regime by subtracting the lensing distortion induced by 
the smooth cluster and galaxy halos. 
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where is the sum of the shear contribution at a 

given position j from A'gai galaxies. This entire inver- 
sion procedure is performed within the lenstool utili- 
ties developed by Kneib (1993), which accurately take into 
account the non-linearities arising in the strong lensing 
regime. Using a well-constrained 'strong lensing' model 
for the inner-regions of the cluster along with the aver- 
aged shear field and assuming a known functional form for 
p{ts) from the field, the likelihood for a guessed model is 



(7) 



We compute C by assigning the median redshift corre- 
sponding to the observed source magnitude for each ar- 
clet. The best fitting model parameters are then obtained 
by maximizing the log-likelihood function with respect to 
the parameters cr* and rt*. Using a bootstrap method, 
we have verified that while the convergence to the best-fit 
model is indeed driven by the brighter cluster galaxies, no 
single galaxy dominates in the procedure. 



3. PROPERTIES OF A 2218 

The lens model for the mass distribution in the well- 
studied cluster A 2218 at z = 0.17 is arguably the best 
constrained case presently available (Kneib ct al. 1996), 
to which we can apply our analysis. We use four sets 
of multiple images, which are all spectroscopically con- 
firmed at z = 0.70, 1.03, 2.52 and 5.60, in conjunction 
with the background shear field, to obtain the best-fit pa- 
rameters for a fiducial galaxy, vt* = 40^^q kpc and 
(T* = 180^20 km s-i (where the quoted error bars are at 
5(7 confidence). The mass distribution in A 2218 is com- 
posed of two clumps centered around the two brightest 
cluster galaxies, and the profile inferred from lensing ob- 
servations and X-ray data are consistent. Although the 
mass distribution is bimodal, we note that the primary 
clump contributes ~ 90% of the mass. Nevertheless, we 
take the presence of the second clump into account when 
calculating the center-of-mass of the cluster. Below, we 
tabulate the model parameters for these two components 
in our best-fit mass model of A 2218. 

The morphologies and luminosities of the cluster mem- 
bers (which are crucial for the lensing analysis) are also 
well-determined from the high-resolution HST imaging 
of this cluster (Small et al. 2001). We have used the 
40 brightest early type-cluster galaxies in the maximum- 
likelihood analysis to obtain the values quoted above. 
Their individual truncation radii were then determined by 
utilizing the scaling laws with luminosity in equation (2). 
Out of the early-type cluster galaxies selected for the lens- 
ing analysis, ^ 75% are located inside the core-radius and 
the rest lie interior to 300 kpc, within a few core radii from 
the cluster center. 

There have been several dynamical studies of A 2218 
that yielded accurate measurements of galaxy velocities 
relative to the center of mass of the cluster (Danesc, De 
Zotti & di Tullio 1980; Le Borgne et al. 1992; Girardi et 
al. 1997; Rakos, Dominis & Steindling 2001). Using the 
spectro-photometric survey of Le Borgne et al. (1992), we 
find that 25 of the 40 early-type galaxies in our analysis 
have measured redshifts which allow us to estimate their 
velocity relative to the center-of-mass of the cluster. We 
have then computed the expected truncation radii for these 
galaxies in two limiting cases regarding the composition of 
their dark matter halos: (i) coUisionless dark matter; and 
(ii) fluid-like dark matter (see Furlanetto & Loeb 2002). 

In our analysis, we have to additionally assume that the 
orbits of the cluster galaxies under consideration are repre- 
sentative of tidally stripped galaxies. Ghigna et al. (1998) 
have found that most of the strongly tidally-stripped halos 
move on nearly radial orbits (since these orbits penetrate 
the high density core of the cluster) and that radial or- 
bits are overabundant by an order of magnitude relative 
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400 ± 30 
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to tangential orbits in the inner r <^ 500 kpc of the cluster. 
Since we are using the projected location of cluster galaxies 
in conjunction with their orbital velocities at the present 
time, the nature of the orbits interior to 300 kpc needs to 
be established in order to justiiy the assumption that these 
galaxies are indeed populating representative orbits for ef- 
ficient tidal stripping. Using the determined mass profile 
of the cluster, the line-of-sight velocity dispersion profile 
and the number density profile of the galaxies, and assum- 
ing that cluster galaxies are good tracers of the potential 
well, wc can compute the velocity anisotropy parameter'^, 
/?(r) = (1 — a^/a1), via the anisotropic Jeans equation us- 
ing the approximation of spherical symmetry (for further 
details see Natarajan & Kneib 1996), 



G Mtot(r)z^, 
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dr 



(8) 



where i'g{r) is three-dimensional galaxy density profile 
(as derived from the 2D projected distribution via the 
Abel transform), (T^(r) the radial velocity dispersion of 
the galaxies, and Miot{r) is the total mass profile. The 
observed line-of-sight velocity dispersion profile ajos {R) is 
given by, 



The above two integro-differential equations can be solved 
numerically for and /3(r). Substituting the best-fit mass 
model for A 2218 we obtain (3 « 0.3, which indicates that 
the majority of orbits are in fact radial and do plunge 
through the core of the cluster. Such orbits are likely to 
facilitate tidal truncation as assumed in our analysis. 

4. RESULTS 

For coUisionless cold dark matter (CDM) the tidal trun- 
cation radius, rf., for the z-th galaxy whose pericenter lies 
at a radial position r,; inside the cluster, is defined by the 
condition that the average density of the galaxy halo inte- 
rior to the tidal radius, matches the mean density of the 
cluster interior to the galaxy location n (Ghigna et al. 
1998; Taylor & Babul 2001): 

{P9Aru)) = {Pci{ri)). (10) 

For the case of fluid- like dark matter (FDM), the halo 
truncation process is modified due to the additional pro- 
cess of ram-pressure stripping. In this case the truncation 
radius is set by the condition (Furlanetto & Loeb 2002), 



(11) 



where Vg^ is the velocity component of the i-th galaxy rel- 
ative to the cluster center-of-mass, PgX^t) is the density 
of that galaxy halo at its 3D truncation radius rj, and 
cjg. is the internal velocity dispersion of that galaxy^. In 
our analysis, we have conservatively used only the mea- 
sured component of the velocity of each galaxy (relative 
to the cluster center) along line-of-sight. Since the full 
(3D) galaxy velocity can only be larger than our adopted 
value, our constraint on the truncation radii of galaxies in 
the FDM case should be regarded as conservative upper 
limits. The considerations quoted above ignore any fur- 
ther truncation that is likely to be caused by thermal or 
turbulent heating (Gnedin & Ostriker 2001). Using the 
above expressions, we derived the truncation radii for the 
25 early-type galaxies in A 2218 in the two limiting cases. 
The results are plotted in Figure 2. The observed distri- 
bution of truncation radii from the Icnsing data analysis 
(solid triangles) is compared to the radii expected in the 
coUisionless limit (solid squares) and the fluid limit (solid 
circles). We find that FDM models are ruled-out at a 
confidence greater than 5a", while coUisionless dark matter 
models are in excellent agreement with the observation- 
ally determined values. We find no correlation between 
the truncation radii and the orbital velocities for the clus- 
ter members in A 2218, as would be expected if ram pres- 
sure stripping had been a significant dynamical process in 
operation. 

The association of the retrieved model radius rt* from 
the lensing analysis with the tidally truncated radius is 
strengthened through the following observed correlation. 
For a galaxy described by a PIEMD mass profile with ve- 
locity dispersion ctq. and whose orbital pericenter crosses 
the cluster core where the density is po, the tidal trunca- 
tion radius scales as 



rt oc pq 



(12) 



Deriving the best-fit values of po from the strong lensing 
models and (Tq. from the maximum-likelihood analysis, we 
find that the above simple scaling adequately reproduces 
the observed trends both for the truncation radius and 
the mass of a fiducial halo in the five HST clusters studied 
in NKS02 (ranging in redshift from z = 0.17-0.58). The 
trends seen in halo size r^ with redshift from the lensing 
study of NKS02 are also in good agreement with the theo- 
retical expectation from numerical simulations (Ghigna et 
al. 1998). High resolution coUisionless N-body simulations 
of cluster formation and evolution, predict that the dom- 
inant interaction is between the global cluster tidal field 
and individual galaxies in the redshift interval z ~ 0-0.5 
(Ghigna et al. 1998; Moore et al. 1996). 



^Note that /3 = implies isotropic orbits, < /3 < 1 implies mostly radial orbits and /3 < indicates orbits that axe primarily tangential. 
^Note that while the truncation radii inferred from the lensing observations are projected 2D values, those calculated from equations (7) 
and (8) are in 3D. However, because the PIEMD mass density profile is p oc {r^ + rQ)~^{r^ + r^)~^, both truncation radii are equivalent. 
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5. CONCLUSIONS AND DISCUSSION 

We have used gravitational lensing data to constrain the 
truncation radii of galactic dark matter halos in the cluster 
A 2218. Figure 2 shows that the inferred truncation radii 
are consistent with the tidal radii expected for collisionless 
dark matter, but rule-out fluid-like dark matter for which 
ram pressure stripping is effective. 

The transition between the collisionless and coUisional 
regimes is set by the ratio between the mean-free-path of 
dark matter particles, A, and the radius, rt, of the galaxy 
halos under consideration. This ratio is given by. 



where ap/ nip is the coUisional cross-section per unit parti- 
cle mass and S(rt) is the surface mass density of a galaxy 
halo at its truncation radius. The fluid regime is obtained 
for — < 1 . The characteristic surface mass densities of the 

rt ^ 

galaxies in A 2218 can be directly inferred from the anal- 
ysis of the lensing data. We find that for an L^, galaxy, 
^*{rt) ~ 0.024 g cm~^. Since the fluid regime is ruled out, 
we exclude all values of (Jp/nip > 42 cm^ g~^. 

Our constraints are complementary to those derived by 



Gnedin & Ostriker (2001) from considerations of ther- 
mal conduction in the mildly coUisional regime. These 
authors exclude the regime 0.3 ^ (Jp/rUp ^ lO^'cm^g"^, 
based on the consideration that cluster ellipticals will oth- 
erwise deviate from the fundamental plane beyond the ob- 
served scatter. Our new constraint allows us to nUe out 
the high cross-section regime and hence we conclude that 
(Tp/nip ^ 0.3cm^g~^ as an absolute upper limit. Dark 
matter cross-sections higher than this upper limit were 
postulated by Spergel & Steinhardt (2000) in order to rec- 
oncile problems that cold dark matter models possess when 
compared to observational data (such as the abundance of 
galactic sub-halos and the slope of the inner mass den- 
sity profile of galaxies; see further discussion in Dave et 
al. 2001; Yoshida et al. 2000; Stoehr et al. 2002; Miralda- 
Escude 2002). 



PN acknowledges support from a Trinity College Re- 
search Fellowship. This work was supported in part by 
NSF grants AST-9900877, AST-0071019 for AL. JPK 
thanks the CNRS and the TMR-Lensing collaboration and 
IRS the Royal Society and the Leverhulme Trust for sup- 
port. 



REFERENCES 



Binney, J., & Tremaine, S., 1987, Galactic Dynamics, (Princeton: 
Princeton U. Press), Ch. 7 

Brainerd, T., Blandford, R., & Smail, I., 1996, ApJ, 466, 623 
Couch, W.J., Barger, A.J., Smail, I., Ellis, R.S., Sharpies, R.M., 1998, 
ApJ, 497, 188 

Dancsc, L., dc Zotti, C, & di TuUio, G. 1980, ApJ, 82, 322 

Dave, R., Spcrgcl, D., Steinhardt, P. J., & Wandclt, B. 2001, ApJ, 

547, 574 

Deir Antonio, I., & Tyson, J. A. 1996, ApJ, 473, L17. 

Fischer, P., ct al. 2000, AJ, 120, 1198 

Furlanetto, S., & Loeb, A. 2002, ApJ, 565, 854 

Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., & Stadel, 

J. 1998, MNRAS, 300, 146 

Girardi, M., ct al. 1997, ApJ, 490, 56 

Gnedin, O., & Ostriker, J. P. 2001, ApJ, 561, 61 

Griffiths, R. E., Casertano, S., Om, M., & Ratnatunga, K. U. 1996, 

MNRAS, 282, 1159 

Kassiola, A., & Kovner, I. 1993, ApJ, 417, 474 

Kneib, J-P., Ellis, R. S., Couch, W., Smail, I. R., & Sharpies, R. 
1996, 471, 643 



Le Borgne, J-F., Pello, R., & Sanahuja, B. 1992, A&AS, 95, 87 
McKay, T. et al. 2002, ApJ, submitted {astro-ph/0108013) 
Miralda-Escude, J. 2002, ApJ, 564, 1019 

Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A. 1996, 

Nature, 379, 613 

Natarajan, P., & Kneib, J-P. 1996, MNRAS, 283, 1031 
Natarajan, P., & Kneib, J-P. 1997, MNRAS, 287, 833 
Natarajan, P., Kneib, J-P., Smail, I., & EUis, R. S. 1998, ApJ, 499, 
600 

Natarajan, P., Kneib, J-P., & Smail, I 2002, ApJ, submitted 
Rakos, K., Dominis, D., & Steindling, S. 2000, A&A, 369,750 
Smail, I., et al. 1997, ApJS, 110, 213 
Smail, I., et al. 2001, MNRAS, 323, 839 

Spergel, D., & Steinhardt, P. J. 2000, Phys. Rev. Lett., 84, 17, 3760 
Stoehr, F., ct al. 2002, preprint, astro-ph/0203342 
Taylor, J. E., & Babul, A. 2001, ApJ, 559, 716 

Yoshida, N., Springel, V., White, S. D. M., & Tormen, G. 2000, ApJ, 
535, L103 

Wilson, G., Kaiser, N., Luppino, G. A., & Cowie, L. L. 2001, ApJ, 
555, 572 



5 




^ Q Q u \ I I \ \ I I I \ \ I I I \ \ \ I I \ \ I I \ \ I I I \ \ I I I \ \ I I I \ \ I I I \ I I \ 

10 20 30 40 50 60 70 80 90 100 

r_t* (kpc) 

Fig. 1. — Results of maximum-likelihood analysis of lensing data for the cluster A 2218 (NKS02). The contour plot shows the best-fit values 
for the model parameters <t» and rt*, which are the central velocity dispersion and truncation radius for a typical L* galaxy in the cluster. 
The likelihood contours are plotted in intervals of 1 — cr starting from the inside out. The thick open curves are lines of constant enclosed 
mass. 
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Fig. 2. — Distribution of truncation radii as inferred from the lensing analysis of 25 galaxies superimposed on the scaling relation (thick 
solid line) in the cluster A 2218. The 3-a and 5-a lines denote the corresponding confidence levels in the parameters rj,, and cr* obtained 
for a typical L, galaxy in this cluster (derived from the confidence level contours in Figure 1). The expected distribution of tidal radii for 
fluid-like dark matter (FDM; solid circles) and coUisionless cold dark matter (CDM; solid squares) are also shown. Note that the linear 
relation between truncation radius and velocity dispersion for the values inferred from lensing is a consequence of the assumed scaling laws 
with galaxy luminosity [equation (2)]. 



